A model for the fragile-to-strong transition in water 
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A model based on the existence of two different competing local structures in water is described. 
It is shown that it can explain the transition between fragile and strong behavior that supercooled 
water has around 220 K. The high temperature behavior is similar to that observed in standard 
fragile glass- formers. The strong behavior at low temperatures is associated with the existence of a 
remanent configurational entropy coming from the possibility of locally choosing between the two 
possible environments. 



I. INTRODUCTION 

The form of the dependence of viscosity on tempera- 
ture is among the many properties that make water an 
anomalous fluid. Water is a fragile fluid when viewed 
at temperatures close to the melting temperature, indi- 
cating that there is an arrest of its degrees of freedom 
on cooling. This behavior is typical of many substance 
known as fragile glass- formers [0J|]. However, close to the 
glass temperature Tq (~ 136 K)supercooled water shows 
characteristics of a strong liquid , in which there is an 
almost temperature independent configurational entropy, 
that manifests in an Arrhcnius dependence of the viscos- 
ity 77 as a function of T. Thermodynamic constraints 
limit the transition between these two regimes to occur 
rather sharply in a temperature range around ~ 220 K 

There is by now a good piece of evidence that many 
of the anomalous properties of water can be rationalized 
by the use of an effective, two particle, spherical inter- 
action potential, of the core-softened type |j] 0|. This 
interaction can be viewed as appearing between clusters 
of water molecules, rather than between single molecules 
pO| , pl| . The main characteristic of this interaction is that 
it allows for two different equilibrium distances between 
clusters, depending on pressure. An appropriate, simpli- 
fied model that capture many of the anomalies of water 
is provided by spherical particles interacting through a 
potential consisting of a hard core plus a soft repulsive 
shoulder |^,^|. Here we show-using an analytically solv- 
able version of it-that this kind of interaction can also 
explain the non-standard behavior of rj(T). 



II. HARD SPHERES MODEL 

We will use a model of hard spheres as a starting point 
(in the next section it will be generalized to describe the 
properties of water). We will suppose that the pure hard 
sphere system has an ideal thermodynamical glass tran- 
sition at some temperature To when the fluid phase is 
supercooled, preventing crystallization. A possible sce- 
nario for this glass transition is the following jl2|-[l4j . For 



glassy systems there is a contribution s c to the entropy- 
referred to as configurational entropy-that comes from 
the many different configurations in which the glass can 
exist. For the case of spheres it comes from the many 
ways in which the spheres can be accommodated in sta- 
ble, non-crystalline arrangements. These configurations 
differ in the value of the specific volume v. We will sup- 
pose that hard spheres have a configurational entropy per 
particle s^ s (v) of the form[| 

s? s (v) = a(v-v )-/3(v-v ) 2 . (1) 

According to this formula, s^f s becomes lower than zero 
for v < vo (and v > vq + a/ft), indicating that there 
are no accessible states in this range, i.e., vq is the mini- 
mum value that v can take. To get the total entropy Stot 
of the system we should still include the contribution 
coming from small vibrations around each configuration 
fL4}| . However, for our purpose this will not be neces- 
sary, since (as we will see below) 77 is given in terms of 
s c , rather than s-tot> and s c (T) is independent of the vi- 
brational contribution to the entropy. The previous form 
of the configurational entropy implies the existence of an 
ideal thermodynamical glass transition occurring at To, 
where T is obtained from P/Tq = dsf s / 'dv\ v _ = a. 
Whereas from microscopic grounds there is no rigorous 
prove that this transition should occur, the consequences 
on observable magnitudes that can be predicted from it 
are consistent with the known phenomenology of glassy 
systems p"5[ and with results of numerical simulations 
Jl^ , ^6[ . For T < Tq the system is in the fundamental 
configurational state, and then s^ s (T < Tq) = 0. For 



T > Tt 



0, s 



is given by 



sf S {T > T ) 



APT 2 



(T2 ~ p2/a2)= 4^ (T2 " T ° 2) 

(2) 

This expression for s^ s can be used to calculate trans- 
port properties such as the viscosity 77 through the use of 



*In Ref. [[[4] a parabolic form for s c as a function of density 
(instead of v) is used. The difference between both choices is 
tiny, and both give rise to an ideal glass transition. 
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the Adam-Gibbs formula |17|. It states that the value of 
rj is given by 



s c (x,v) = {l-x)s? b (v) + 



V (T)=r, exp[A/(Ts c )}, 



(3) 



which for hard spheres becomes 

V HS (T) = 770 exp [4Aa- 2 (3T/(T 2 - T 2 )] (4) 

where r/Q and A are constants. This is a behavior typ- 
ical of a fragile glass-former, in which r; increases more 
rapidly than in a simple thermally activated process, and 
it diverges when T — » To. The presence of s c in (^|) re- 
flects the fact that jumps between different basins of the 
energy landscape become less probable as the number of 
these basins diminishes. 

III. CORE-SOFTENED MODELS FOR WATER 

Properties of water have been studied recently by us- 
ing models in which particles interact through potentials 
that allow for two different equilibrium distances between 
particles, namely do an d di > do ||-§,[l8). One possibil- 
ity is, for instance, to take a strict hard core at do and 
a shoulder that vanishes at d\ HQ. Here we will use a 
further simplification of this kind of models in order to be 
able to extract analytical results . We consider spheres 
of radius r\ (= d\/2"). Pairs of spheres will be allowed 
to overlap (more than two overlapping spheres will not 
be allowed), and each time this happens the system will 
be charged an energy £o. This may be considered as a 
limiting case of particles with a core-softened potential, 
in which there is a low energy hard core at a distance 
2ri, and a strict hard core at 2ro, and we are taking 
ro = 0. To make the problem analytically tractable, we 
will also suppose that each time two spheres overlap, they 
are constrained to have their centers in exactly the same 
position. This approximation neglects the entropy asso- 
ciated with small vibrations of the particles in each pair. 

We are interested in the configurational entropy s c of 
the system, now as a function of the specific enthalpy 
h = Pv + e, that includes the internal energy e coming 
from the existence of overlapping particles. To calculate 
s c (h) we proceed in the following way. Suppose we have 
a system of iV particles, n of them in non overlapped 
positions and n' pairs of overlapped particles (N = n + 
2n'). The configurational entropy will be that of n + n' 
hard spheres plus the combinatorial entropy for locating 
the n' pairs in the n + n' possible positions, i.e., 



n HS 



N 



ks In 



n + n 
n' 



(5) 



(here we use s c to indicate an entropy functional). Using 
x = n'/N and v = V/(n + n') = v/(l—x) as independent 
variables we can write s r as a 



(1 - 2x) In 



1-x 
l-2x 



■ In 



(G) 



In order to get the thermodynamic value of s c (h) at each 
value of the external pressure P, we have to maximize s c 
for each fixed value of the enthalpy, i.e., 



s c {h) 



mm 

x,v 



s c (x,v) 



h = Pv + e = (1 — x)Pv + x£q. 



(7) 
(8) 



IV. RESULTS 



In Fig. H we see s c , x, and v as functions of h for three 
different values of P, using expression (^) for s HS (v) with 
a = 2.79 k B r^ 3 , = 0.97 k B r^ 6 , and v = 6.37 r\ which 
are values extracted from numerical simulations of hard 
spheres systems Q. We also plot in Fig. |l] the limiting 
cases s c (x = 0) and s c (x — 0.5), corresponding to all 
particles in singled or overlapped positions, respectively. 
The thermodynamic value s c (h) can never be lower than 
s c (x — 0) or s c (x = 0.5). The states with the lowest 
enthalpy for x = and x = 0.5 have (from (ih) h = Pvq 
and h = Pvq/2 + Eq/2, respectively. These values coin- 
cide at P cr = eq/v . In Fig. |(a), P = 0.9 P cr , S c (x = 0) 
is always greater than s c (x = 0.5) and for this reason 
x takes values close to 0, indicating that most particles 
are in singled positions. The obtained s c (h) function de- 
parts from zero with infinite derivative at h = 0.9 £o (be- 
cause of the combinatorial contribution to the entropy) 
but it still has a singularity (namely a jump in its sec- 
ond derivative) when ds c /dh = a/P. In (b) the value of 
P = 1.1 P cr is larger, and at low h the contributions with 
x = 0.5 dominate, indicating that the system has almost 
all particles coupled in pairs. For higher h, x goes down 
to zero, namely paired particles become rare. In (c), the 
value of P — Per is exactly that at which the ground 
state of the system with x — and the one with x = 0.5 
are degenerated. In this case, entropy starts from a finite 
value ~ 0.48 ^ at h min — Sq, corresponding to the max- 
imum combinatorial entropy of choosing which particles 
are singled, an which are paired. 

To get the values of the thermodynamic variables as a 
function of T, instead of h, we have only to make use of 
the relation = ds/dh. The results for the configu- 
rational entropy s c (T) and the viscosity r](T) (calculated 
using the Adam-Gibbs formula (|J)) are shown in Fig. ||. 



V. DISCUSSION AND COMPARISON WITH 
WATER 

From Fig. |^(a) we see that for any P, s c is finite for 
all T / 0, i.e., there is no vanishing of s c at any finite 
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temperature, contrary to what happened in the case of 
simple hard spheres (see Eq. (||)). This is due to the 
possibility for the system of having particles singled or 
paired, which always accounts for the existence of a non- 
zero combinatorial entropy. For P — P CT the bottoms 
of the s c functions corresponding to x = and x = 0.5 
coincide (Fig. 0(c)), and this combinatorial entropy can 
be used up to T = 0, in such a way that s c remains finite, 
even when T — > 0. For P 7^ P cr , s c goes to zero as T — > 
0, as the ground state is unique. There is still a phase 
transition at finite temperatures, signaled by the kink in 
the s c (T) curves in Fig |^(a). The position of this kink in 
the P-T plane is given by P/T = a. In the log(ry)-l/T 
plot (Fig. ||b)), the kinks mark the transition between 
a fragile and a strong behavior. At high T, s c diminishes 
rapidly with temperature and the system is fragile. This 
behavior is equivalent to that of the simple hard sphere 
system (see Eq. (Q)). At low T, the dependence of s c on 
T is much weaker, indicating a stronger behavior. In this 
regime, the available configurational entropy is mainly of 
combinatorial nature. 

For this model there is no ideal glass transition, i.e., 
r\ is finite at any finite temperature, and Tq = 0. How- 
ever, from an experimental point of view, the glass tran- 
sition Tq is conventionally defined as the value of T at 
which 77 takes some large value (this is usually taken to 
be 10 13 poise). In Fig. || we show curves of constant r\ 
extracted from our model. Each of these may be thought 
as defining a dynamical glass transition temperature Tq 
(that depends on pressure), according to different dy- 
namical criteria. The position of the line of the fragile- 
to-strong transition is also indicated. We see that Tq is 
systematically lower around the critical pressure P = P CI . 
This behavior has been observed in numerical simulations 
of SiC>2 [ ^9|j20| , which has a fragile-to-strong transition 
qualitatively similar to that of water. Notice that Tq is 
lower than the temperature of the fragile-to-strong tran- 
sition only in some range of pressure around P cr , and 
only in this range the fragile-to-strong transition will be 
experimentally observable. 

One important ingredient that we have not included 
in the model is the existence of an attractive part in the 
interaction potential. In real water this attraction gen- 
erates the liquid-vapor first order coexistence line, and 
also probably a second first order line in the supercooled 
region separating two different amorphous configurations 
. A simple way of analyzing the consequences on our 
model of an attraction between particles is the following. 
If the attraction is considered to be long ranged, of van 
der Waals type, then all the results we have obtained re- 
main valid if we replace P by a new effective pressure 
P* = P + 7/w 2 with some constant 7, namely, the at- 
traction acts as an effective pressure (which depends on 
v) that has to be added to the external pressure P. In 
the P-T phase diagram, this non-uniform transformation 
of the P axis produces (if v decreases rapidly when P in- 



creases) a "folding" that indicates a first order transition 
. This is the way in which the liquid- vapor coexistence 
line appears in the van der Waals equation for a fluid. 
For our model, since v(P) at T = has an abrupt dis- 
continuity at P = P cr , the attraction generates also a 
new first order line ending in a critical point. It seems 
to be pj that in water this critical point lies at temper- 
atures higher than Tq, so the first order line determines 
two different sectors of the glassy phase of water. These 
two sectors correspond to the experimentally observed 
high density amorphous (above the first-order line) and 
low density amorphous (below the first order line) phases 
|Io| , pT| . In our model, these two phases differ in the frac- 
tion of particles that are paired, and thus they can also be 
named high density and low density amorphous phases. 

We have relied for our discussion upon the existence of 
a thermodynamic phase transition for hard spheres, that 
is not rigorously proven to occur. However, for slightly 
different forms for sf? s (which may imply the absence 
of an ideal glass transition) than that given by Eq. (||), 
our results still remain valid, except for the fact that the 
sharp fragile-to-strong transition (the kinks in the curves 
of Fig. ||) becomes a crossover. 

The model we have presented explains the transition 
between fragile and strong behavior of water as appearing 
from the competition between two different local struc- 
tures. To be able to solve the problem analytically, we 
had to make the crude approximation that these struc- 
tures correspond to singled and paired particles. In real 
water it is likely that what plays the role of our parti- 
cles are the so-called Walrafen pentamers , which are 
clusters of five water molecules. These cluster are nat- 
urally expected to accommodate at one of two possible 
distances from each other . It is clear that in this more 
general case the mechanism for a fragile-to-strong tran- 
sition may remain basically the same. In fact, from the 
way we solved the model, it is seen that all what is needed 
is the existence of two different competing structures, in- 
dependently of the details of them. At high T the fragile 
behavior is associated to the configurational entropy of 
each structure individually. At low T the strong behav- 
ior appears due to the combinatorial entropy of choosing 
locally between the two structures. 
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FIG. 1. Configurational entropy s c , specific volume v, and 
fraction of pairs of particles relative to the total number of 
particles x as a function of the enthalpy for three different 
values of pressure. In order to compare, the limiting cases 
s c (x = 0) and s c (x = 0.5) are also shown (entropies are given 
in units of ks, v is in units of 5rf, where n is the radius of 
the particles, see text for details). 
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FIG. 2. (a) Configurational entropy s c as a function of T, 
for different values of P/P CT , as indicated, (b) Viscosity r\ vs. 
1/T, calculated according to the Adam-Gibbs formula (ex- 
pression (Eh, r/o and A are the constants in that expression). 




FIG. 3. Pressure-temperature plot, indicating the position 
of the fragile-to-strong transition (solid line) and contour lines 
of constant viscosity (dashed lines, the corresponding values 
of EoA -1 \n(rj/rjo) are indicated). 
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